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Abstract 

Self - consistent temperature dependence of the average magneti- 
zation in quantum Heisenberg ferromagnet is obtained as a first ap- 
proximation of perturbation theory on an inverse radius by apphcation 
of the functional method for quantum ferromagnet. The experimental 
data are compared with the theoretical results for EuO ferromagnet. 
A quantitative agreement is observed between theoretical and exper- 
imental results outside interval near Curie temperature. The theory 
is of critical behavior in the vicinity of Curie temperature Tc- The- 
oretical Tc value itself is agreed with the experimental one with an 
accuracy of several percents. 

1 Introduction 



The functional representation of a partition function allows to extend the 
application of the approximate methods to the problem of the statistical 
physics. Convenience and naturalness of the functional representation re- 
sulted in numerous attempts to extend this method to the Heisenberg quan- 
tum ferromagnet model (see 0, ^, Q). A method proposed in |]^ makes it 
possible to write an expression for the partition function of a magnet in the 
asymmetric phase in the explicit closed form of an integral over two number 
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fields — a neutral field and a charged field. The charged field \1/ corresponds 
to the elementary excitations of magnet. The neutral field p corresponds to 
the average spin magnitude. Erroneous interpretation of global effects in 
was corrected in 0]. 

The partition function may be considered as a generating functional of 
averages of the spin operators. Accordingly the closed expression for the par- 
tition function allows to calculate approximately different averages of spin 
operators, using small parameter 01/2; — the inverse range of the interac- 
tion. 

A problem of the average spin calculation can be formulated within the 
scope of the functional formalism as a minimum of the effective potential 
Weff{p) problem |^. We will calculate Weff{p) in the first order in 1/z. 
Minimizing We//(p), we will next obtain the average magnetization M{T) as 
a function of the temperature. Such a way of considering it allows to avoid 
the emergence of the fictitious divergences presented in the direct diagram 
expansion. The comparison between theoretical relationship M(T) and ex- 
perimental one for the ferromagnet EuO demonstrates that they are in close 
agreement, just as magnitudes of Curie's temperature with an accuracy of a 
few per cents. 

2 Generating functional of averages of the 
spin operators 

The model of the Heisenberg quantum ferromagnet is based upon the Hamil- 
tonian 

H = --^SiJijSj, (1) 

where Jij — is the matrix of the exchange interaction between neighboring 

spins Si which are located at the lattice sites r^. If the interaction of spins 
with external magnetic field is taken into account the partition function can 
be written as: 

Z[hi{t)] = TrTexp | -/3H + J dthi(t)Si | , 
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where the symbol T denotes time ordering and hi{t) is the external field at 
lattice site fi. For the operator exp{—eH) a Gaussian transformation 

expi-eH) = I n^^(detJ,,)-i/2gxp \-Y,$^J^'^,e/2 + $Sie\+0{e^) 



can be performed to within terms ~ e^. Thus, writing exp(—PH) = (exp(—eH))^^'^ 
with e — > 0, we arrive at an expression for the partition function in the fol- 
lowing form 
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X 



HTrTexp [ | dt [0;(t) + hi(t)]li 



(2) 
(3) 



(4) 



The substitution put forward by Kolokolov in P| 

(f>! = JijPj - hi 

<Pt = i^t - hf, 

rearranges the ordered operator exponential (|]) to the explicitly specified 
operator 



A{(3) =Texp 



dt [<P{t) + h{t)]S = exp(V^-(/3)5+) 
X exp ( S' f dt pit)] exp f exp ( f p{t') dt'] dt 



(5) 



xexp(-V^~(0)^+), 



where index i is omitted, 
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Thus, regarding ip^{t) and p{t) as new integration variables, we can calculate 
the trace of the T-exponential @ explicitly and obtain a closed functional 
representation for Z[h]. We have for the arbitrary spin magnitude S 



2S 
1=0 



\p-{f3)-^'{0)) f dt^+(t)exp / p(t')dt' 



/ /3 



p(t) dt 



m=~-S m=—S 

Pi = S{S+l)-{m + l){m + l-l). 

The authors of showed that it is necessary to use the initial conditions 
V'j~(0) = in order to transform from the original variables to the new 
integration variables (^). With this conditions and in the discretization of 
the transformation fffl) 



lim - \ ip~{t + e) — ip' 



(6) 



dt e 

the Jacobian of the transformation is equal to the constant. Hence follows a 
functional representation for Z[h] 




(7) 



\o 



'13 



gi \rp^{t)dt \ I (8) 



+ 



-hi) J, 



(9) 



This partition function may be considered as a generating functional of av- 
erages of the spin operators: 

5hi(z\h]) 
< Stit) >= lim =< P^it) > ■ 



shm 
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3 Effective potential of field p 



In order to obtain the field p effective potential , we must integrate over 
fields ip^. We cannot perform this integration exactly, however it can be 
made approximately over parameter 1/z = J2j Jfj / JijY = Sfc "^1 / "^o; 
where is the Fourier transform of the exchange matrix Jjj = J(rj — r^) 

(jo = j(^ = o), EfcJfc = o). 

Let us expand the nonpolynomial part of action r„ @ in a series in 
powers of ip'ip^ and extract all the F terms at /i = 0, which are linear 
functions of 



Fl = / dtY.\pr{t)Jr,p,{t) + 2^+(t) ^ - i,-{t)J,,p,{t) + 

ij \ / 

/ dtY.2b{pm) (^+(t)V^rW -C(/3)V^+Wexp(pO(t))no(p°(/3))) - 







llgMm. (10) 

i 

t 

P°i{t)=IEdT JijPj (r) , no (0 = exp(i)-i ' 

j 



where = (iln(5fo(0)/'^^ Brillouin function. The contributions from the 
other terms of F are small in the inverse range of interaction, because these 
terms include two and more sums of {JijY- Hence in a zero approximation 
on 1/2; follows: 

1. The saddle-point value p is determined by equation 

<^^>^/ = 6(/3Jo/), (11) 
which agree with mean field theory equation. 

2. The Fourier transform of the longitudinal correlator takes the form 



(12) 
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where p = j^J2Pi = p{k = 0) is the average value of the field pi, M — 

i 

number of lattice sites. 

3. The bare propagator of the field ip^ in the Fourier representation has 
the form 



Gdt, t') =< i^^imU^') >^ y (^(^ - + + 1) hp(j£&t) - 1]) X 

(13) 

exp(^^(t-t')), 

n-k = r~m 7' = ^oP - Jkb, b = b{pJop), 9{0) = 1, 

exp(a;^/5) - 1 " 

where 6{t) — theta function. 

4. The transverse correlator in the Fourier representation has the following 
form 

= 4Jr2 < > -2J^'5it - t') ^ 

FX/ FX/ FX/ FX/ 

2b (n,^ +l-9it- t')) exp(^g(t - t')). (14) 

Let us remind that we impose the zero initial conditions ip~{0) = on 
the ip~it) field, then the propagator is of the G^(0,t') = property and 
is nonperiodical. Function tug does not approach zero value at = for an 
arbitrary value of p. However, if p is defined by (|TT]) in the zero approximation 
on 1/z, the Goldstone symmetry will be restored (u;(0) = 0). 

By substituting for F„ + F; and integrating (|^) over ip"^, we obtain the 
partition function in the first approximation on 1/z. 

Z = N fUDp, exp{-W\p,)), (15) 

where up to a constant W^{pi) has the form F^ in which the bare propagator 
(|13D is substituted for ip~ (t)ip~^ {t') . It can be pointed out that W^{pi) also 
determines a longitudinal correlator, whose infrared behavior was examined 
in low-temperature limit in 0. 
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Let us perform the Fourier transformation of field pi and expand into 
a series in p{k 7^ 0). It is enough to hold the quadratic terms in p[k 7^ 0) 
only within the first order in 1/z. By integrating over p{k 7^ 0), we obtain 
the effective potential W^ij)). The saddle-point value p is determined in the 
first approximation on by minimizing the effective potential W^ij)): 

P = & - E(l - PWn^ - no) + ^ E T^^' (16) 

k k ^ 

This equation was obtained in by spin diagram technique application, 
wherein (|T6|) was concerned as the approximation to (|Tl|) . Within the scope 
of our approach the equation (|T^) means the equation in the saddle-point 
value p in the effective potential, so there is a need to solve ([l6|) exactly for 
calculation < >= p. The temperature dependence of the S^(T) at S'^ — > 
can be found analytically. Magnetization S{T) approaches zero as a square 



S{T) = A {tI - t) ^ (17) 



at temperature 




k 



(18) 



T° = 6'(0) Jo = ^^^Jo, 
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which is the critical temperature value in the first approximation on 1/z. 

For the illustration of the results obtained, let us compare magnetization- 
temperature relationship ([16| ) with that obtained in experiment. We have 
pick Heisenberg ferromagnet EuO as a standard for comparison. Its struc- 
ture (FCC — face-centered cubic lattice), exchange integral values for inter- 
action between neighboring spins Jl = 1.212''i^' (12 neighbors) and J2 = 
0.238°i^' (6 neighbors), magnetization-temperature relationship S{T) and 
Curie's temperature = 69.15^ K were examined carefully in experiments 
[|]. In EuO the magnitude of spin is equal to 5* = 7/2 and the perturbation 
theory parameter is HkJ^ / Jl = 1/14, Jq = 15.97°/^. The figure illus- 
trates experimental relation between temperature and magnetization S{T) 
(points) and relation given by mean field theory equation (0) (curve 2) and 
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given by equation (|T6j) (curve 3) which were calculated by computer. Similar 
comparison between theoretical results obtained by spin diagram technique 
and experimental data on EuO and EuS were made in ||^. The jump of 
curve 3 at 65^ K determines critical region. The method of steepest descent 
with the help of which the equation ([T6| ) was obtained, does not work in the 
critical region. It means that dependence S(T) specified by ([1^) breaks down 
over critical region in the vicinity of T^. However, Tc itself can be obtained 
exactly, since the effective potential dependence W{p) changes qualitatively 
at T = Tc. For T < Tc function W{p) has two minimums at p = ±p, and for 
T > Tc function W{p) has one minimum at p = 0. Hence, calculating effec- 
tive potential at each following order of perturbation theory on 1/z, we can 
obtain more exact equation in magnetization as well as critical temperature 
value. 

It is significant that the critical behavior which is characteristic to the 
phase transitions of the second kind in the method of effective potential is 
conserved in the perturbation theory. Therefore, in spite of inapplicability of 
our approach near Tc, it nevertheless does not disturb the qualitative picture 
of the phenomenon. 

The author is thankful to I.V. Kolokolov for numerous discussion and a 
help, and also to V.S. L'vov for interest to the work. 
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